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Abstract 

Principles of discrete time mechanics are applied to the quantisation of Maxwell's 
equations. Following an analysis of temporal node and link variables, we review the 
classical discrete time equations in the Coulomb and Lorentz gauges and conclude that 
electro-magneto duality does not occur in pure discrete time electromagnetism. We dis- 
cuss the role of boundary conditions in our mechanics and how temporal discretisation 
should influence very early universe dynamics. Quantisation of the Maxwell potentials 
is approached via the discrete time Schwinger action principle and the Faddeev-Popov 
path integral. We demonstrate complete agreement in the case of the Coulomb gauge, 
obtaining the vacuum functional and the discrete time field commutators in that gauge. 
Finally, we use the Faddeev-Popov method to construct the discrete time analogues of 
the photon propagator in the Landau and Feynman gauges, which casts light on the 
break with relativity and possible discrete time analogues of the metric tensor. 

Throughout this paper the acronym CT refers to continuous time, whereas DT refers to 
discrete time. Readers familiar with the principles and methodology discussed in the earlier 
papers of this series may skip the introduction, but are advised that the general notation 
has been improved and is discussed in section 2. 



1 Introduction 

The development of quantum field theory in the second quarter of the twentieth century was 
accompanied by speculation about the microscopic nature of time and space, motivated by 
the spectacle of quantised fields dynamically evolving over a classical space-time in no way 
different to the bland Riemannian space-time continuum used by Einstein in general relativ- 
ity. Various problems in quantum field theory such as the divergences in the renormalisation 
programme and ambiguities in operator products were believed to be associated in some way 
with the microscopic description of space-time, but relatively little was done to investigate 
the issue in any depth. Although there were occasional attempts to investigate alternative 
mathematical descriptions of space-time, such as the notable work of Snyder fl], §, it was 
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more usual to circumvent difficulties with ad hoc procedures such as point splitting of oper- 
ator products with no modification of the underlying space-time, or by the use of space-time 
lattices, which were always assumed to be an approximation to the continuum. It was only 
with the advent of the space-time foam approach to quantum gravity and the reinterpre- 
tation of superstring theory in the nineteen eighties that it became generally acceptable to 
talk about the microscopic nature of space-time as more than likely quite different to the 
continuum normally assumed in field theory. It is in this context that our work should be 
seen. 

In the current series of paper on DT mechanics [3-6] we investigate the consequences of 
taking literally the hypothesis that time is discrete on an incredibly small scale. Our original 
motivation is discussed in the first paper of this series |§. Compared with quantum gravity 
and superstring theory, ours is a very modest and limited step. However, it turns out to have 
enormous consequences, principally conceptual, altering virtually all aspects of the laws of 
mechanics and how we view space-time. Along the way a number of sacred cows have to be 
sacrificed. For instance, without continuity with respect to time there is no differentiation 
with respect to time. Therefore we are forced into the construction of a mechanics without 
velocities. This immediately raises the question of what replaces Lagrangians, which are 
normally functions of dynamical variables and their temporal derivatives. Since we have 
no velocities in our theory, we cannot construct conjugate momenta in the traditional way, 
as these are defined as derivatives of Lagrangians with respect to velocities. So we appear 
not to have a phase space, and consequently we do not have a Hamiltonian formulation or 
Poisson brackets in the normal sense of the word. This goes hand in hand with the lack of 
continuous translations in time and with the absence of a generator of such transformations. 
This raises questions about quantisation, but we have shown in earlier papers of this series 
that these can answered. The price we pay for this is that we end up no longer doing exactly 
what we were doing before, amounting to a paradigm shift in the language of Thomas Kuhn. 

We regard it as an important principle that we do not simply modify CT equations of 
motion by replacing temporal derivatives with ad hoc differences. Although that works in 
some situations in Newtonian mechanics, it becomes more subtle in the presence of gauge 
invariance. Our approach is to start from the beginning, rewriting the action integral as an 
action sum and developing the consequences rigorously from there. In this series of papers 
we have not confined our interest to particular models which happen to be amenable to 
temporal discretisation. The microscopic nature of space and time will affect all dynamical 
processes. Our interest has been in the subject of DT mechanics as a whole, in both its 
classical and quantised forms, applied to point and field systems. 

An important consequence of temporal discretisation is that it rejects the notion that 
space and time form a four dimensional continuum. This notion has been successfully ex- 
ploited in the special and general theories of relativity throughout this century and is a 
cornerstone of the theories of quantum gravity and superstring theory (which paradoxically 
eventually undermine this very idea). We appear to have taken a step back towards the 
separation of absolute space and time in Newtonian mechanics. Discretising time also raises 
questions any relativist would ask, which are: in which inertial frame is time discrete and 
what dictates this choice? Our discrete time mechanics is not Lorentz covariant, and we 
have to address this issue as well as others. 

We believe that there are satisfactory answers to these particular questions which accord 



2 



with modern cosmology. Consider a CT Freidmann- Robertson- Walker space-time, i.e., one 
for which CT co-ordinates can be chosen so that the Riemannian pseudo-metric distance 
rule takes the standard form 

ds 2 = dt 2 -R(tfda 2 . (1) 

Here R (t) is a function of the co-moving time t only and da is the distance element of 
some spherical, flat, or hyperbolic 3— space. Given that the gross space-time structure of the 
expanding universe is well represented by such a choice of co-ordinates and by such a distance 
rule, consider now the cosmic microwave background radiation (CMB) field discovered by 
Penzias and Wilson in 1964 0. It has been pointed out by a number of authors ||, [J, 
that, contrary to the principles of relativity, this radiation field can be used to define a 
local absolute inertial frame Fp at each point P in space-time. Such a frame is unique up 
to Euclidean transformations such as rotations. For an observer at P instantaneously at 
rest relative to Fp, the CMB radiation field will appear isotropic to a very high degree. 
That this is physically meaningful is borne out by the empirical observation that the earth 
appears to be moving at a speed of about 500 — 600 km/sec relative the local Fp frame, a 
phenomenon called the dipole effect |10| . It is somewhat ironic that long after the Michelson- 



Morley experiment failed to detect an aether carrying radiation (and thereby supporting the 
principle of special relativity), Penzias and Wilson discovered a plenum consisting entirely 
of radiation which may be used to define absolute local inertial frames. 

Our thinking is as follows. Suppose we take seriously the hypothesis that time is really 
discrete. Our fundamental criterion is that any dynamics based on this idea should not make 
predictions at odds with scientifically determined (i.e. empirical) facts. We shall call this 
the empirical principle. It has been part of our programme to determine where if anywhere 
temporal discretisation actually clashes with this principle. So far, we have not been able to 
rule DT mechanics out on this basis. 

Always mindful of this principle, we should now carefully sift out and identify those 
additional concepts and ideas which, despite our traditional training and inclination, are 
really no more than contemporary belief structures which happen to be compatible with the 
empirical principle. We shall call such ideas idealisations. We should feel free to dispense with 
idealisations if absolutely necessary, provided we do not clash with the empirical principle. 
It is particularly important to identify idealisations which are currently popular because of 
their elegance and mathematical content; these are usually the hardest to remove, because 
mathematical elegance is frequently taken as a principle in physics. 

A good example of such an idealisation is the Poincare group. This particular structure 
cannot occur in our theory because we do not have the freedom to make continuous dis- 
placements in time; this property emerges only in the continuous time limit. What replaces 
the Poincare group in mechanics is still under investigation. What can be said how- 
ever is that any DT analogue will almost certainly be more complicated mathematically. 
The absence of the Poincare group however does not unduly worry us here and should not 
be used as a criticism, because our mechanics can still satisfy the empirical principle for 
the following reason. In our DT mechanics there is a fundamental interval of time T and 
generally, we can show agreement with CT theory certainly at the O (T°) level. Moreover, 
disagreement invariably occurs at the O (T 2 ) level. Since we imagine T is of the order of 
the Planck time T P = ^/hG/c 5 ~ 5 x 10~ 44 sec, our theory should be good for the current 
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level of experimental accuracy. As yet we have no explanation of the value of T, but then 
we cannot explain the value of c, H, G or the electric charge either. 

Continuing this line of thought, we argue that temporal discreteness would have influ- 
enced both the dynamical origin of the universe and its subsequent evolution. The Planck 
epoch is a term used to denote the interval from the origin of time to Tp. By the end of that 
interval, it is generally believed that gravity had decoupled from the other interactions and 
space started to expand in a pre-inflationary context. It is during the Planck epoch that 
CT field theories are generally considered to be either invalid or seriously incomplete, and 
it is possible that some version of DT mechanics (though not necessarily the one we are 
considering here) is the appropriate theory to describe very early universe dynamics. If true, 
then a generally covariant description in the fashion of general relativity might not be at all 
appropriate during this epoch, with the best description perhaps involving some preferred 
frame of reference. 

It may be the case that after the Planck time and before false vacuum inflation, different 
regions of the Universe had different local temporal discretisation frames, randomly dis- 
tributed, analogous to pre-inflationary early universe monopoles or ferromagnetic domains 
in a solid. If so, we would argue that in much the same way as the monopole problem was 
removed by inflation, only one actual temporal discretisation frame would survive inflation 
into our local universe. All the others would be beyond the event horizon. 

The discretisation frame holding in our visible universe would most likely be linked to the 
dynamical processes involved with the phase transition from the false vacuum to radiation 
and matter, and we would expect it to leave a signature. This signature would be the 
existence of the co-moving frame used in the FRW metric (|1|), and subsequently, the isotropy 
frame of the CMB. We propose that the discrete time parameter discussed in this paper be 
identified as the discretised version of the local co-moving time (assumed to be the coordinate 
time in our local Fp) in our neighbourhood. We are assuming here that the CMB isotropy 
frame and the cosmic matter rest frame coincide locally |5|. 

Now over laboratory scales associated with particle scattering experiments, we could 
ignore the local aspect of this discretisation and regard a global discrete time frame as a very 
good approximation, in much the same way that special relativity is a good approximation 
to general relativity in the laboratory. Hence we end up with the equivalent of a unique 
discretisation of Minkowski space-time in our neighbourhood. 

A DT analogue of general relativity awaits investigation; it will require abandoning the 
principle of general covariance on a microscopic level, with a generally covariant description 
emerging only in the CT limit. It is possible that our fundamental interval T is itself 
determined dynamically by the local matter densities and by the dynamics. There is a 
precedent for this idea. In 1983 Lee published a study of discrete time field theory [|TT|j , 
where his interval of time T n varied and was a dynamical variable coupled to matter. In our 
earlier work and in this paper, time is discrete but otherwise is as passive as time in CT 
relativistic field theory. We work with a fixed T throughout. 

In this paper we develop further the discrete time Maxwell's equations first discussed in 
|3|, [|]. In the next section we review our notation, which has been somewhat overhauled 
and compactified compared to earlier papers in the series, but with no change in content. 
Following that, we review the DT mechanics formalism, extending the discussion to include 
the dynamical variables associated with temporal links. Our view of discrete time mechanics 
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has evolved from our original picture of dynamical variables changing over successive instants 
of time (nodes) only to one where there are dynamical variables on the nodes and on the 
intervals or links between the nodes. A gauge theory such as Maxwell's equations (and 
undoubtedly gravitation) involves a dynamical interplay between node variables and link 
variables. 

Then we review the DT Maxwell's equations in the Lorentz and Coulomb gauges. Work- 
ing in the former gauge presents a particularly interesting challenge in mechanics because 
it is a relic of the Lorentz symmetry which DT mechanics undermines. At first sight is ap- 
pears less natural than the Coulomb gauge defined in the local absolute rest frame in which 
time is discretised. We recall that the use of Lorentz gauges and attempts to maintain 
manifest Lorentz covariance in CT quantum electrodynamics leads to problems with state 
vectors, requiring the Gupta-Bleuler formulation or equivalent technology. 

We then consider the quantisation of the free Maxwell fields. Using the DT Schwinger 
action principle which was successfully used previously for the DT quantisation of the scalar 
and Dirac fields |6j , we find the vacuum functional and the free field commutators in the 
Coulomb gauge. Although we do not have a Hamiltonian formulation in our theory, we do 
have gauge symmetry, and we may use the Faddeev-Popov approach to the quantisation of 
gauge fields in our theory. We find that the Faddeev-Popov vacuum functional for the DT 
free electromagnetic theory in the DT Coulomb gauge coincides precisely with the vacuum 
functional obtained by solving the DT Schwinger functional differential equations in the 
same gauge. 

Finally, and this was found to be the hardest task, we construct the DT analogues of the 
photon propagator in the Landau and Feynman gauges, in preparation for future applications 
to QED. We find that there is an interesting breakdown of what happens in CT theory, with 
DT analogues of the metric tensor making appearances. These differ in interesting ways, 
depending on whether we are in space time or in momentum space, and indicate that a 
reformulation of general relativity via DT mechanics will be quite instructive. 



2 Notation and conventions 

In addition to the natural unit system where c = H — 1, we use the following conventions 
throughout this paper. Given T is the fundamental time step in DT mechanics then an 
event in our space-time has natural coordinates x = (n, x), where n is an integer. The CT 
limit 

T -> 0, n -> oo, nT = t (t fixed) (2) 

then corresponds to the CT coordinate time t. Our coordinates are natural in the following 
sense. We work in an inertial frame of reference such that the CMB dipole anisotropy is 
absent locally |I(J. In this frame we fix cartesian spatial axes with coordinates x = (x 1 , x 2 , x 3 ) 
such that the distant stars do not appear to rotate relative to such axes. Finally we choose 
an instant of discrete time as our temporal origin of coordinates and count fundamental 
intervals of T forwards and backwards, giving us the integer coordinate n referred to above. 
Our discussion in this paper does not include any cosmological aspects such as Hubble 
expansion, other than using the CMB to determine our frame of reference. We shall be 
concerned with local coordinates relevant to particle physics in the present epoch. 
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We denote integration over space-time by the sum/integral symbol 



rf 3 x. (3) 



Given a variable f x = f n (x) indexed by a discrete temporal index n and by a continuous 
spatial index x we define the DT fourier transform f p of f x by 



f P = f (z, P) = f e px f x , (4) 

J X 

where 

e px = z n e~^ (5) 

and z is complex and non-zero. We shall assume that such a transform exists and defines a 
function analytic in some annular region in the complex- z plane centred on the origin and 
including the unit circle. This of course puts restrictions on the sequence defined by f n . 
Although p is the direct equivalent of the spatial components of the momentum space four- 
vector of CT mechanics, the analogue of the time component p° of the momentum space 
four vector is not z. We shall be interested in taking z on the unit circle, and then it is the 
principal argument of z which is related to p°. We shall use the symbol p to denote the set 
p = (z,p) . 

If the DT fourier transform exists and has a Laurent expansion in some annular region 
centred on the origin in the complex z— plane and containing the unit circle, then we can 
construct the inverse transform. This is given by 



fx = f e px f p , (6) 
J P 

where 

with the z— integral being over the unit circle in the complex plane taken in the anticlockwise 
sense and 



n ip-x 



^px 

We could attempt to define our dynamics in the (z, p) transform space, restricting our 
dynamical variables to be functions in the complex z— plane for which the inverse transform 
(§) exists. No new dynamical content should emerge from this approach, but it may be 
mathematically more secure. However, we will normally assume that coming first from the 
space-time direction is valid. That after all is the conventional way to define field theories. 

The analogues of the identity operators (delta functions) in our formalism are defined as 
follows. If 5 n is the DT Kronecker delta, satisfying 

5 n = 1, n = 0, 

= 0, n^O (9) 
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where n is an integer, then we define the four dimensional DT Dirac delta 5 X by 

^ = %^ 3 (x-y), T >0, (10) 

where y = (m, y) . Then 

$ fx$x-y = fy (11) 

J X 

An integral representation of 5 x - y is 

$x—y Y ^-px&py (12) 
Jp 

The corresponding operator 5 p _ q in DT fourier transform space satisfies the relation 



fpSp-q = f q , g=(«,q), (13) 

p 

where u is complex. K DT sum/integral representation of 5 p - q is given by 

Sp—q X Cpx^qx- (1^) 

The classical step operator U n acting on any temporally indexed variable f n is defined 

by 

U n f n = /n+1, /n = /n-lj (15) 

with powers of U defined in the obvious way, viz 

U a J n = f n+a: (16) 

where a is an integer. We note that 

U n fm = fm, m ^ 71. (17) 

Throughout our mechanics we shall deal with operators and equations defined via linear 
combinations of the step operators. An operator of the form 

P (U n ) = c UZ + Cl U« +l + ... + c r U^ r (18) 

will be called an r th order operator. Then an equation involving an r th order operator will be 
called r th order. Important first order operators are the forwards and backwards differences 
defined by 

A: ee U n - 1, A" = 1 - U' 1 (19) 
respectively, and from them we define the second order symmetric difference 

A n ee \ (A+ + A;) = \ (U n - U- 1 ) . (20) 
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First temporal derivatives are invariably replaced by one of three possible operators, defined 
by 

D+ = ^n. d~ = ^ D = — (21) 
whereas the second derivative is invariably replaced by the operator 

Dl ee DID- = D- n Dt = Un - 2 T + U "\ (22) 

It is a feature of our DT mechanics that the formalism will tell us which of the above 
difference operators we need to use in a given context. For example, our investigations into 
the DT Schrodinger equation || and the Dirac equation show that in those equations 
we have to replace d/dt by the second order symmetric operator D n . This has important 
consequences as far as the solutions of the equations are concerned, because the second 
order symmetric difference leads to a DT equation of motion which acts like a second order 
CT equation of motion, rather than a first order equation of motion, and this generates 
the oscillon solutions discussed in 0] and ||. Fortunately, we found that in the second 
quantised theory, these oscillons correspond to states with unphysical norm, and so are not 
observable asymptotically as ordinary particles. Their role as virtual particles in QED will 
be considered in the next paper of this series. 

In our mechanics, time reversal amounts to the interchange U n <-> U^ 1 . We shall en- 
counter various difference operators which have the symmetry property that they are invari- 
ant to time reversal. We shall refer to such operators as T— symmetric. They are important 
and useful to us, and in fourier transform space they are real functions of z provided 
z is on the unit circle (we note that z* = z~ l holds only on the unit circle). With this 
definition, we see D 2 n is a second order T-symmetric operator. Another important second 
order T-symmetric operator which occurs throughout our mechanics is given by 

5 n = i(C/ n + 4 + [/- 1 ). (23) 

The factor of 1/6 and the 4 can readily be understood in terms of our virtual path procedure, 
discussed below. In the CT limit (fj) , if it exists, these operators can be replaced by 

D+ -> d u D-^dt, D n ^dt, (24) 
Dl - dl (25) 

There is no trace of S n in CT theory, but it plagues DT mechanics, occurring in unpre- 
dictable places and introducing temporal nonlocality in unexpected places. It makes the 
rewriting of CT mechanics into a DT form far from easy, particularly in the case of Maxwell's 
equations. This nonlocality also enters at the level of the metric tensor, and suggests that 
an attempt to formulate general relativity into a DT framework will require thinking of 
the metric tensor as a non-local in time operator (and hence as a more dynamical object), 
rather than as a set of local functions forming the components of a rather bland second rank 
tensor. Another important operator is the DT d'Alembertian, which turns out to be the 
second order T-symmetric difference-differential operator 

n x ee D 2 n - S n Vl , (26) 
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which appeared in our studies of the DT Klein Gordon equation |4], [5|] and the DT Dirac 
equation || . This operator is also important in our formulation of Maxwell's equations. 
We note the useful result 



f x P(U n )g x =j f x P(U- l )g x . (27) 
From this we see that for T-symmetric operators, we may write 

f,P(U n )g x = j f*P(U n )g x . (28) 
The DT fourier transform gives the following useful result: 

e px P (U n ) f x = P (z- 1 ) f p , z ^ 0. (29) 



We may use the above to find the DT Feynman propagator A Fx , which satisfies the 
equation || 

(D x + m 2 S n ) A Fx = -5 X . (30) 
Taking the DT fourier transform of this equation we find 

(p 2 - m 2 Sz) Ap p = 1, (31) 

where we define 

p 2 = -D 2 z - Szp.p, (32) 

with 

z-2 + z- 1 z + 4 + z- 1 
Dz= — , Sz = . (33) 

A solution of interest in particle theory is 



A Fx = <b e px - 2 , (34) 

J p p z — m z bz + le 

where we choose a DT analogue of the Feynman +ie prescription. The singularity structure 
in the complex z-plane of the integrand in the above is particularly interesting. First, we 
can prove that the equation 

p 2 — m 2 Sz + ie = (35) 

has no solution on the unit circle in the complex z— plane for any value of the linear momen- 
tum p. To prove this, write z = exp (i6) . Then fl3"5l) becomes 

2 (cos^ -1) (cos# + 2), 2 v . n 

-±-^r 2 L - '- (p • P+m 2 ) + ie = 0, (36) 

which has no solution for real 9 if e > 0. This result is important because it means we have 
a fully closed contour of integration over the unit circle in the complex z— plane, requiring 
no principal value discussion. 
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Next, we note that (^) is a quadratic in z, with roots Z\, z 2 satisfying the relation 

Zl z 2 = 1. (37) 

Hence we deduce that the denominator in ( |3~3| ) contributes one simple pole inside the unit 
circle and one outside. By looking carefully at the location of these poles as the spatial mo- 
mentum p varies, we see that two distinct patterns of behaviour emerge |3j . For momentum 
in the elliptic regime, corresponding to momenta bounded by T|p| < vT2, the simple pole 
inside the unit circle is just inside and gives the equivalent of a trigonometric solution when 
we use the calculus of residues. For momentum in the hyperbolic regime, on the other hand, 
given by T|p| > vT2, the simple pole interior to the unit circle starts to move towards z = 0, 
giving a damped exponential solution when we use the calculus of residues. 
Finally, we note that if we had taken a slightly different prescription, viz 

Afx = / e — — — — (38) 

J P p — yrn 1 — xe) Sz 

it is not hard to see that our conclusions would be exactly the same. 



3 Link and node variables in DT mechanics 

The discretisation of time requires us to change the way we think about dynamical variables. 
Given a lattice structure to time, we can identify two distinct geometrical components. These 
are the node sites, which correspond to events at times nT, where n is an integer, and the 
links between these nodes. For a one dimensional lattice, links and nodes are mathematically 
dual, but in our theory this mathematical duality does not carry over into the physics. Some 
of our dynamical variables are defined at nodes whereas others are defined at sites, and 
generally this occurs in such a way that interchanging links and nodes is not a symmetry 
of the theory. In particular, matter fields corresponding to massive particles are defined on 
nodes only. We have found that electric fields occur as link variables whereas magnetic fields 
are node variables, so that electro-magnetic duality does not appear to occur here. This 
seems a natural way to explain the absence of Maxwellian magnetic monopoles. 

We show now that the dynamical rules for link variables are similar to but not identical 
to those for node variables. First, let a generic node field dynamical variable be denoted 
by A® (x), where n is the time, x is the spatial position, and a is some extra label such as 
a vector, spinorial or group index. Likewise, denote a typical link variable by </>^(x). As 
discussed in P, || we base our dynamics on the system function F n = J d~x.J :n , where the 
system junction density T n has the form 

<f> n , VA, VA n+1 , V(j> n ) , (39) 

where we have suppressed the spatial coordinates and field labels. This is the discrete time 
analogue of a Lagrangian of the standard form 

L = dx£ (<p, V<^, <p) (40) 
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and is equally generic, in that all of the system functions we need to use are of this form. 
We note that link and node variables are treated differently right at this point, in that the 
system function is first order in the node variables but zeroth order in the link variables. 
This emphasises further that as far as dynamics is concerned, link and node variables are 
not dual. 

In DT mechanics the action integral becomes an action sum. The action sum from initial 
time MT to final time NT > MT is given by 

N-l 

A NM = T ^ p n ( 41) 

n=M 



and the equations of motion are obtained by Cadzow's action principle fT2| 

5A NM = 0, (42) 

c 

for suitable variations of the fields. Here and elsewhere we shall use the symbol = to denote 

c 

an equality holding by virtue of the equations of motion. For an arbitrary variation: 

A° - AZ + 6AZ, M <n < N 

^ n - 0j + <ty£ M<n<N (43) 

we find 

saNM = ? fd 3 x{ 5A a M (x) -J—-F M + y 6A« (x) -^--F n 
J \ 5A a M (x) n ^ +i ; M-(x) 

+M; ' (x »sw f " + £ M » (X) ^ F " _I (44) 

iv \ i n=M+l n v ' 

N-l r 




n=M ^ n v ' 

where we use the summation convention for the field labels a only. For fixed end-point, but 
otherwise arbitrary variations, viz. 

6A a M (x) = 6A% (x) = (45) 

we find the functional derivative equations of motion 

= . , r \F n + F n ~ 1 } = 0, M<n<N (46) 

d "" 0, M<n<n } (47) 



c 



which reduce to 



{^n + ^n- 1 } = V- — — — -{r n + ^ n - 1 }, M < n < N (46 



M«(x) J c <9VA-(x 



= V ■ RVAfif ^ M < n < N. (49) 

Z (x) c 5V< (x) 
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We see here once more the essential difference between the node field dynamics and 
that of the link fields. The former involve genuine second order difference equations of 
motion ( f48| ) corresponding to second order temporal derivatives in CT theory. For the links 
however, equations (|49D are at most first order difference equations, and such equations are 
analogous to constraint equations in CT mechanics rather than equations of motion. All of 
this is intimately tied in with the boundary conditions required to solve these equations. CT 
theories with constraints and gauge symmetries will have DT analogues where both sort of 
equations occur. The electromagnetic equations discussed in this paper give a basic example 
of these ideas. 

4 Classical electromagnetism 
4.1 The charge free Maxwell equations 

Our discrete time formulation of the charge free Maxwell's equations starts with the CT 
electromagnetic potentials A^ 1 = (0, A) which are used to construct the physical electric and 
magnetic fields E and B. A clear distinction has to be made here between the nature of 
the electric scalar potential and the magnetic vector potential A in DT mechanics. The 
former is associated with the temporal interval or link connecting times t n = nT and t n +i, 
whereas the latter is associated with nodes, or times t n themselves. This distinction also 
manifests itself in the difference between the physical electric field E and the magnetic field 
B, which are likewise associated with temporal links and nodes respectively. The electric 
scalar potential associated with the link connecting time t n and t n+ i at spatial position x 
will be denoted by the symbol <p n (x) , rather than by (say) <p n+ i (x) . Although our notation 
suggests a bias towards t n at the expense of t n+ i, this is not really the case. 

One of the principles we have applied to DT mechanics is minimality. By this we mean 
our belief that if we are really working at an irreducibly fundamental level, then the dynamics 
should be as basic and uncomplicated as possible. Now as discussed in a system function 
is more like a Hamilton's principle function than a Lagrangian, being a function of the 
dynamical degrees of freedom at the end points of a fundamental interval of time T. In 
principle, we cannot probe below this level in DT mechanics. We could simply postulate a 
system function, but a better approach would be to take some CT theory such as Maxwell's 
electromagnetism and construct a suitable system function directly from its CT Lagrangian. 
If we did not use the CT Lagrangian as a guide, then we could not hope to guess an 
appropriate system function. 

Our approach uses the auxiliary concept of virtual path. This is discussed fully in [||, |[] . 
The basic idea is to replace the CT fields by functions of the node and link variables smeared 
in some way over a given interval [nT, (n + 1) T] . From the point of view of CT mechanics, 
this introduces a degree of non-locality in time, but this is inevitable. For some fields, such 
as the neutral scalar field, the function is a linear interpolation, but in the case of matter 
fields involving gauge invariance, such as Dirac fields, the virtual paths are highly non- 
linear. Fortunately, the virtual paths for the Maxwell potential fields are straightforward. 
The virtual path nA for the electric potential (a link variable) is defined by 

nA (x) ee n (x) , 0<A<1, (50) 
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where we use the variable A to interpolate between the two ends of a given link. Given that 
the magnetic vector potential is a node variable, its associated virtual path A n \ is defined 
by 

A nA (x) = AA n+1 (x) + AA n (x) , 0<A<1, (51) 

where A„, A n+ i are the dynamically meaningful field values at the ends of a given link and 
A = l-A.. 

The above virtual paths are defined over the interval [nT, (n + 1)T] . For a set of successive 
time intervals, we observe that virtual paths for node variables are continuous in time, 
whereas virtual paths for link variables are not. This is tied in with the nature of their 
respective dynamics. Note that no physical meaning can be attributed to the virtual path in 
a genuine DT theory, and as we have stressed before 0, these paths are used only as a tool 
in obtaining a system function with built in properties such as gauge invariance. Beyond 
that, they have no significance. 

The DT version of Maxwell's equations comes with a DT version of gauge invariance. 
A local DT gauge transformation involves a gauge function x associated with nodes rather 
than links. A gauge function value at time n and position x will be denoted by x n ( x ) an d is 
assumed differentiable with respect to x. The virtual path for a gauge function is given by 

X nA (x)EEA Xn+1 (x)+A Xn (x), 0<A<1. (52) 

In CT electromagnetism a local gauge transformation is defined by the replacements 

A» _> A' 1 * = A» + d^x- (53) 

In our theory this becomes 

0n -» (P'n = K + D nXn, (54) 

A n -> A; = A n -Vx n . (55) 

When we talk about local gauge transformations, we shall mean these last two equations. 

Turning now to the physical fields, we define the gauge invariant electric and magnetic 
fields via the potentials: 

E n EE -V0 n - D+A n , B n = V x A n . (56) 

These satisfy the homogeneous DT Maxwell equations 

V • B n = 0, V x E n + D+B n = 0. (57) 

To construct a gauge invariant system function, we recall that the CT Lagrange density 
for electromagnetism is given in terms of the Faraday tensor F^ v = d^A v — d u A lx , viz 

C — —-F 

= -\d^A v & i A v + \d v A^A v 

= liA' + d^iA' + d^ + ^d^d^-d^diA"). (58) 
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The gauge invariant system function density for the charge free system is obtained by replac- 
ing the fields in the above CT Lagrange density by their virtual path forms and integrating 
with respect to A over the interval [0, 1]; 

T n = (C(A nX ,<P nX )) 

= \ (D+A n + V0 n ) • (D+A n + V0 n ) + '-((dtAi&A^ - d.A^A^)), (59) 
where the angular brackets denote an integral over A, i.e. 

(/ A >= / f(X)d\. (60) 
Jo 

Here we use the virtual path replacement dt — > T~ 1 d\. For convenience we have not multiplied 
the system function by a factor T as was done in earlier papers of this series, so that it now 
has the physical dimensions of a Lagrangian rather than an action. 

The equations of motion for the magnetic potential (a node variable) are given by 

d - t {T n + r 1 - 1 } = d 3 -^- {t- + r*- 1 } , (6i) 



dA i J c J djdA 

which reduce to 

□ n A n + VA„ = 0, (62) 

c 

where 

A n = D-<P n + V • S n A n (63) 

is the DT Lorentz function. For the scalar potential (a link variable), the equation of motion 
is 

v ' im: (64) 



c 



which reduces to 



This implies the equation 



D+V-A n + V 2 0„ = 0. (65) 

c 



□ n n - D+A n = 0. (66) 

c 



Equations (p2| — 166|) are DT gauge invariant, as can be readily verified. The reason is that 
our system function is gauge invariant, and this guarantees that the equations of motion are 
gauge invariant. 



4.2 The DT Lorentz gauge 

In the DT Lorentz gauge we set 



An = 0, (67) 
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and then the potentials satisfy the massless DT Klein-Gordon equations 



□ n A n = 0, (68) 

c 

□ n 0„ = 0. (69) 

c 

It is always possible to work in this gauge, as the following argument shows. Suppose 
that we started off with a configuration of fields for which the DT Lorentz function was 
non-zero, i.e. 

A n = D~(j) n + S n V-A n ^ (70) 
Now consider the gauge transformation 

A' n = A n - VX„, <P' n = K + D nXn (71) 

where the gauge function x n satisfies the DT inhomogeneous Klein-Gordon equation 

U nXn = -An. (72) 

Then we find 

A' rt = D~4f n + S n V-A' n = 0, (73) 

as required. We note that using our experience with the DT Klein-Gordon propagators in 
earlier papers, we may write down a particular solution to ( [72]) in the form 



X, 



f^F.-yK (74) 
J v 



where A Fx is the DT scalar massless Feynman propagator satisfying the equation 

n x A Fx = -8 X . (75) 

It is noteworthy that the Lorentz condition 

A n ee D-<f> n + S n V-A n = (76) 

is second order in time, that is, relates field values on three successive nodes and on the two 
links between them. This will interact in some way with the equations of motion, which are 
also second order, and so we can expect trouble. We shall see that the Lorentz condition and 
its generalisation demands special attention when we come to work out the electromagnetic 
propagator in the DT Feynman gauge. In particular, we will have to consider the existence 
of an inverse operator, S"" 1 , which is highly non-local in time. 



4.3 The physical fields 

Turning to the physical (gauge invariant) electric and magnetic fields, we may write the 
system function density ( [BHD in the form 

jF n = ^E 2 n — - (B n+1 • B n+1 + B n+1 • B n + B n -B n ) , (77) 
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which differs from the more familiar form 

F' n = \^ 2 n — |B n ■ SnB n (78) 

by a total temporal difference and so gives the same equations of motion. Applying Cadzow's 
equation to (ff7|) we find 

V • E n = 0, D~E n = V x S n B n . (79) 

c c 

Using (|57D and ( ff9|) we find the physical electromagnetic fields satisfy the DT massless 
Klein-Gordon equations 

□ n E n = 0, □ n B n = 0. (80) 

c c 

We may readily construct the analogues of the conserved total linear momentum and 
angular momentum using the method described in Appendix B. For example, we find the 
DT analogue of the Poynting vector is 

P n = J d 3 x {E„, x B n + \TB\ V B\ +1 } = P n+1 . (81) 

The expression for the total angular momentum is left as a exercise. 

4.4 Maxwell's equations in the presence of charges 

In the presence of electric charges the system function density ( |77| ) is replaced by 

jrn =T n_ + l_ An+i . + 1 ^ . jn (g2) 

where p n (x) and j n (x) are the discrete time charge density and charge current respectively. 
The homogeneous equations (|57|) remain unaltered but now the equations of motion become 



V • E„ = Pn , V x S n B n - D~E n = j„. (83) 

c c 

These equations are consistent provided the equation of continuity 

^nPn + V-j„ = (84) 

c 

for electric charge holds. The dynamical equations of motion ([83]) may be written in the 
form 



- W„ - D+V-A n = p n . (85) 

c 

□ n A n + VA n = j„. (86) 



Equation ( p5|) may also be rewritten in the form 

□ n n - D+K = S nPn . (87) 
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In the Lorentz gauge (|86|) and ( 



87|) become 



= S nPx , (88) 

c 

□xA x =j x . (89) 

c 

We note here the appearance of the non-local T-symmetric operator S n . 
Quantisation is more convenient in the Coulomb gauge, where we set 

V-A n = 0. (90) 

Then the equations of motion (^) and fl8"6| ) become 

V 2 n = - Pn , (91) 

c 

□ n A n + VD~<j) n = j n . (92) 

c 

We see once again that the scalar potential <p n cannot be regarded as a dynamical field in 
the same way as the components of the vector potential are. Equation ( |9~T| ) is zeroth order in 
time whereas (|92|) as a full second order dynamical equation. This is the DT analogue of the 
situation in CT electromagnetism, where a direct application of Dirac's constraint analysis 
shows that the momentum conjugate to the scalar potential vanishes. This has important 
consequences when we develop our quantisation via the Schwinger action principle, discussed 
next. 



4.5 comment 

It comes as a surprise to see that charge density as formulated in our mechanics turns out to 
be a link variable. A naive guess would have us take it to be a node variable, on the grounds 
that electric charge is carried by matter fields, which in our theory are node variables. This 
is another example where the unravelling of dynamics from a CT to a DT framework forces 
us to re-evaluate our understanding of the various components of dynamics. 

For example, consider the preparation of a state containing charged particles. In view of 
the above comment about charge being a link variable, we see that it must be insufficient in 
some way to think of such a state as being completely defined or specified at a given instant 
or node of time only. If we wanted to measure the total charge of the system, for instance, 
we would have to consider the fields on the link to which this node is attached, and the fields 
on the node at the other end of this link as well. 

This raises another interesting thought; the link pointing forwards in time from a given 
node is different to the link pointing backwards in time. Therefore, the meaning of what 
constitutes a state in DT mechanics must depend on whether it is regarded as an initial 
state or as a final state. 

It is such examples which lead to the conclusion that in cosmology, time could not 
be considered to have a beginning at a point only; it would be necessary to say something 
about the first link as well. If indeed it is correct to think of our fundamental interval T as 
equivalent to the Planck scale Tp, then the origin of the universe in DT mechanics would 
require us to regard the so-called Planck epoch as just the first link. Then it would be wrong 
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in this context to imagine any form of dynamical evolution process occurring during that 
epoch. This is in direct contrast with fundamental theories such as quantum gravity and 
superstring theory based on continuous time, where presumably, there is scope for a great 
deal of dynamical interaction during the Planck epoch. 

Carrying on this line of thought, we would have to accept that dynamical fields on nodes 
such as the vector potentials A could only start dynamical evolution after the Planck epoch, 
since they satisfy a second order equation of motion. On the other hand, link variables such 
as the scalar potential (ft would not have such a restriction. The essential point here is that 
DT mechanics alters our perception of boundary conditions. Fortunately, our discussion in 
this paper involves the present epoch, and we may assume time runs from remote past to 
remote future without any qualms about boundary conditions at the origin of time. 



5 Quantisation in the Coulomb gauge 

Quantisation in field theory is readily tackled via the DT Schwinger action principle, 
which we shall now state and use to determine the DT electromagnetic field commutators 
and vacuum functional in the Coulomb gauge. It is a merit of Schwinger's approach and also 
of Feynman's path integral approach to quantisation that the emphasis is on the physically 
useful amplitudes of the theory, rather than on the operators themselves, such as occurs in 
the canonical quantisation process. The process of imposing naive canonical commutators 
between dynamical variables and their conjugate momenta can be expected to fail in DT 
mechanics for a number of reasons: we do not have a Hamiltonian framework in our mechan- 
ics; we do not have a constraint theory in the fashion of Dirac for gauge field dynamics; and 
the construction of conjugate momenta is straightforward only in the case of systems which 
are normal 0. Fortunately, none of these reasons prevent us from quantising Maxwell's 
equations. 



5.1 The DT Schwinger action principle 

When we use the DT Schwinger action principle, we work in the Heisenberg picture and 
consider matrix elements between states associated with different times. These correspond 
to preparation and observation, that is, initial and final states. If |\&, M) is the state we have 
prepared at time MT, and |$, N) is a state we are asking questions about at time NT > MT, 
then the infinitesimal change N\^>, M) in the transition amplitude ($,iV|^/,M) due to 
infinitesimal changes in the external sources is defined by 

<5($,JV|tf,M) = i{$,N\8A NM \%,M), N > M (93) 

where 5A NM is the infinitesimal change in the action sum operator 

N-i . 

a nm = T J2 d 3 *F n - ( 94 ) 

n=M J 

In the case of the electromagnetic field, we take the free field system function density ( [77| ) 
and introduce arbitrary infinitesimal sources p n , j n in the manner of Schwinger ||13|| . Since 
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these sources are arbitrary, we must take care to ensure that the charges which couple to 
the electromagnetic fields satisfy the equation of continuity ([84]) . Following Schwinger and 
anticipating the use of the Coulomb gauge, the system function density in the presence of 
the sources is given by 

T n [j] =T n - KPn + \ An+1 " 'f n+ l + \K ■ j° , (95) 

where 

Pn( X ) = PnW, (96) 

j°(x) = j n (x) + vjG Cx _ y [D-p y + V y .j y ] (97) 

Jy 

are the conserved charge densities constructed out of the independent external densities p n 
and j n , and Gc is the Coulomb Green's function which satisfies the equation 

V 2 x G Cx = -8 X . (98) 

This has particular solution 

G Cx = (99) 

° 4tt|x|T v ; 

We note the conserved charge densities satisfy the DT equation of continuity 

£>-pP(x)+V x .j«(x) = 0, (100) 

regardless of the values of the independent densities. The above coupling ensures that we 
are free to vary p n and j n arbitrarily whilst still ensuring that the electromagnetic fields are 
coupled to conserved charges. 

Functional differentiation is defined via 

Jify s J^T) P ™ (y) = ^ s (x " y) (101) 



and similarly for the currents. Then with the external sources coupled as in fl95|) we find 

—^—(^,n + l\^,n) j = ($,n + l|0 n (x)|¥,n)j (102) 

5p n ( x ) 

^—(^n^n-l^ = 1 ^,n\Ai^)\^,n-l) j (103) 

-Z—^n+ll^n^ = i($,n + l|4(x)|*,n>,-, (104) 

using the Coulomb gauge (or transversality) condition 

($|V-A n (x)|tf)=0 (105) 

for all states $. 



19 



If now we assume the existence of in and out vacua, and if Z [j] = (0 out \0 in )j denotes the 
vacuum functional in the presence of the external sources, then direct application of the DT 
Schwinger action principle gives the functional derivatives 



id 



-id 
&3n (X 



■Z\j) = (0 out \(f) n (x) |0 in )j, 
■Z[j] = (0^14 (x) \0 in )j. 



In the Coulomb gauge the quantum analogues of Cadzow's equations are 



x 



-Pn W 



a n (Q^\A n (x) |0-), + V xJ D"0 n (x) Z [j] = r n (x) Z [j] 



(106) 
(107) 

(108) 
(109) 



taking the scalar potential in the Coulomb gauge to be a c-number. From these equations 
we deduce 

(110) 



f Gcx-yPy 
J y 



and 



□.{(T'Kitr),- = %z [ 3 ] + d?d*t G Cx . y3 ^z [ 3 ] . 



Note that we have switched notation, as we shall do frequently, using the symbols x, y to 
denote (n, x) and (m, y) respectively. 

The second functional derivative of this last equation gives the DT time-ordered product 



I I \ i n 

(0|fXi4|0) = iA Fx _ y 5ij + ij& Fx _ z d*d*G. 
in the absence of the sources. This is equivalent to 



Cz-y 



<0|2X (x)^ m (y)|0) 



where (p) is the fourier transform 



d 3 p 



ptpj 



A „-m 



P) e 



-ip-(x-y) 



:ii2i 



(113) 



A£(p) 



d 3 xA n F (x) e~ ip - 



(114) 



of the temporally indexed Greens' function, which satisfies the DT massless field equation 

5, 



□ n A" F (x) = -^ 3 (x). 



(115) 



In the absence of external sources we may choose the radiation gauge, defined by 

n (x) = O, ($|V.An(x)|tt)=0, Vtt,$. (116) 
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Using the above time ordered products ( |113| ), we can extract the following unequal-time 
commutators: 



<0| I; +1 (p),4 + (q) |0) 



— iV r . 



(0| AZ, (p) , Ai (q) |0> = -iT 



p 2 

pipj 



(2vr) 3 5 3 (p-q) 
(2tt)V (p-q) 



where 



and 



6T 



4(P) = 
^ + (P) = 

Hence we obtain the result 

(0|[4 +1 (x),i4i(y)] |0) = -i 



6 + T 2 p 2 



^xe-tp-x^ ( x ) 



rf 3 xe ip x ^ (x) . 



(117) 
(118) 

(119) 

(120) 
(121) 



rf3 P -ip-(x-y) 



(27T) 



6T 



6 + T 2 p 2 



:i22i 



which is precisely the same as for the scalar field discussed in |J apart from the modified 
Kronecker delta, necessary to preserve the transversality condition ( |105|) . 

As a final step, we may suppose that the commutators of the fields are c-numbers, in the 
language of Dirac, and then we arrive at the operator commutator statement 



[Ai +1 (x) , Ai (y)] = -i J 



d 3 p 

(27T) 3 



3 ip-(x-y) 



6T 
6 + T 2 p 2 ' 



(123) 



which amounts to our DT quantisation prescription. We may use 
equation of motion 



0. 



and the operator 
(124) 

(125) 
(126) 



a n (p, A) = — J ^xe^-^ e (p,A) • [A„ +1 (x) - e i6 *A n (x)] (127) 
a+ (p, A) ee j (fxe^" e- (p, A) [A n+1 (x) - e -**A„ (x)] (128) 



□ n A n (x) = 0, 



to deduce the equal time commutators 

~K (p) , Ai+ (q) 

which is equivalent to 



[4 (x) , ,4£ (y)] =0. 
Photon creation and annihilation operators are defined by 
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where 

„ 6-2T 2 p 2 , N 

cos9 ^YTtW (129) 

and 

e (p, A) -e (p, A') = 5 XX >, e (p, A) p =0. (130) 

In these definitions we have taken the momentum p to be in the elliptic regime. From our 
work in previous papers with the DT harmonic oscillator, we know that momentum in the 
hyperbolic region |p|T > y/l2 would give wave-functions which either decay to zero or diverge 
in time. The implications are that, like scalar and Dirac particles discussed in earlier papers 
of this series, there is a natural cutoff in our mechanics in the photon spectrum. This may 
have important repercussions in discussions involving for example the black body spectrum. 
A direct application of the commutation rules (|123| — 1126|) gives 



[a n (p, A) , a+ (q, A')] = 2|p|^l - ^jfhx' (2vr) 3 5 3 (p - q) (131) 
[a n (p, A) , a n (q, A')] = 0, (132) 

which shows explicitly that we will obtain a spectrum of polarised photon states, but only 
up to the parabolic barrier |p|T <\/l2, as discussed above. We note that an expansion in 
powers of T gives 

[an (p, A) , a+ (q, A')] = 2|p|5 Ay (2vr) 3 5 3 (p - q) + O (T 2 ) , (133) 

which supports the results discussed in earlier papers of this series that Lorentz symmetry 
emerges very rapidly from our mechanics if, as we imagine, T is of the order of the Planck 
time or less. 



5.2 Comments 

i) Fixing the value of 6 to equal 9 P in ( |129| ) is the DT equivalent to a mass-shell constraint 
in CT field theory; 

ii) The cutoff in the high momentum photon particle spectrum occurs for a dynamical 
reason. Consider the following toy model; suppose ip n is a complex valued variable evolving 
according to the DT harmonic oscillator equation || 

ip n+1 = 2r]ip n - ip n _ x (134) 

c 

where rj is real. Now take another complex sequence (z n ) satisfying the same equation. Then 
we can readily prove that the constructions 

< = -« M*+i-^+iC) (135) 
are invariants of the motion, viz 
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regardless of the value of 77. In a field theory, the equivalent of 77 is determined by the 
momentum p. 

Now we ask where the cutoff in comes from. The answer is found in the pattern of possible 
behaviour found in the DT harmonic oscillator. We showed in that bounded motion 
occurs when 77 is restricted to the range — 1 < 77 < 1, corresponding to the trigonometric 
solutions of the CT harmonic oscillator. We call this regime the elliptic region. For 77 = d=l 
we have the parabolic regime, and for \q\ > 1 we have the hyperbolic regime, where the 
solutions either diverge or collapse to zero in the limit of infinite time. It is this which 
creates the cutoff. If we attempted to create particle states with momentum in the hyperbolic 
regime, we would find physically unacceptable behaviour occurring in matrix elements after 
a sufficiently long time. Such states would not be stationary in the conventional sense. An 
analogous phenomenon occurs in CT quantum wave mechanics, where we reject solutions to 
wave equations on the basis that they have unacceptable behaviour at large spatial distances. 



6 The vacuum functional 

A fundamental problem in quantum gauge field theory is the construction of the vacuum 
functional Z [j] . In this section we shall first solve the functional differential equations 
obtained via the DT Schwinger action principle to obtain the pure electromagnetic vacuum 
functional in the presence of external sources in the Coulomb gauge. A functional integral 
approach then becomes of interest as an independent means of calculating the same quantity 
and confirming that the DT Schwinger action principle is sound. 



6.1 The DT Schwinger vacuum functional 



In the Coulomb gauge we take n to be a classical field and the A n to be quantum fields, 
satisfying the DT Schwinger action principle relations 



iS 



-i5 



If now we take the equations of motion 



% Gcx-yPy 
J V 



and write the vacuum functional as the product 

z [j] = z[ P ]z [j] 



(137) 

(138) 
(139) 

(140) 



then we may readily integrate the resulting functional differential equations to find the DT 
Schwinger action functional 



exp 



xJy 



Px^Cx-yPy + fx^-Fx-yjj 



-i 



Fx — yGrC'y—z 



2.3 



6.2 The DT Faddeev-Popov vacuum functional 



It is a remarkable feature of quantum gauge field dynamics that we may use two completely 
different routes to determine the CT vacuum functional Z \j\ expressed as a path integral in 
CT field theory. Superficially these seem quite different. One way is to work out the details 
of the physical phase space dynamics via Dirac's constraint mechanics and then define Z [j] 
in terms of the non-redundant or physical degrees of freedom. The other route is via the 
Faddeev-Popov gauge symmetry approach. These apparently unrelated approaches give the 
same result. 

The former approach involves the Hamiltonian [14], [15| and therefore has no analogue in 
DT mechanics as far as we understand at this time (because we do not have a Hamiltonian 
approach in our mechanics, there being no generator of continuous translations in discrete 
time). 

Fortunately, the Faddeev-Popov method approach uses gauge symmetry arguments to 
achieve the same end, and there is nothing in those arguments which prevents us from em- 
ploying exactly the same method in gauge theory. By following the standard arguments 
14, 15] we arrive at the expression 

Z [j] ~ II / t rfA «] ^] A 9 [A,0] 5 [g (A,0)] exp {iS (142) 

where g is the gauge fixing function. Here we have taken the gauge fixing function to reside 
on the nodes. This is consistent with the Coulomb gauge discussed next and with the DT 
analogues of the Landau and Feynman gauges discussed after that. 
If we choose the Coulomb gauge then 

9n = V ■ A n , (143) 

the A g factor is independent of the fields, and the action integral reduces to 

S [j] = j {40V 2 - |A • DA - p<j> + j • A} . (144) 

Then using the functional integral result 



J [d-^] exp | z J dx ipMip + i j dxjip^ ~ Vdet M _1 exp | — |i J dxjM ^ j- 



(145) 



we recover the Schwinger vacuum functional ( |141| ) exactly. This confirms the consistency of 
our approach. 



7 More general gauges 

The Faddeev-Popov approach allows us to consider more general gauges. We now follow the 
standard gauge fixing approach to construct the DT analogues of the Landau and Feynman 
gauge propagators. First we note that the vacuum functional ( |142|) is independent of the 
gauge fixing function, so we can choose 

g x = A x - uj x , (146) 



24 



where u x = uj n (x) is an arbitrary function on DT space-time and A is the DT Lorentz 



function (|63|) . Now we can functionally integrate over to, weighting the integrand by a 
suitable weighting function. The essential step here is to choose the weighting function to 
be the exponential 



ox]) | — ^ u x N(U n )u x 



(147) 



where N (U n ) is some T— symmetric operator to be determined and a is a constant. After 
integrating out u via the functional delta function, we can take the vacuum functional to be 
given by 

Z [j] ~ J [dA] [d(j>] exp {iS + iS G + iSj} , (148) 

where So, Sq and Sj are the free electromagnetic action sum, the gauge fixing term and the 
source term respectively. We now consider these three terms separately. The most efficient 
approach is to work with the four component object 



41" 



-Vu A n ) , 



(149) 



which looks like a Lorentz four-vector, but it should be stressed that this is a matter of 
convenience only. For one thing, A° n = (p n is a link variable whereas the other components 
are node variables. 

Suppressing the space-time indices (n, x) , we find 



So = 



^{-f ^DA* + - 2 <),.\ : Si) r \' - |0V 2 + d^D+A*} 



where the operator 



-D-d, 



-OSij — Sdidj 



written in (1, 3) x (1, 3) block form acts to the right. 
The gauge fixing term is given by 



s 



G 



~ j, (D~(j> + Sd.A') N (D~(j> + SdiA 1 ) 



+ 1 -a- 1 ^S^A" 



where the operator 



S, 



ND 2 ND+Sdi 
'" " NSD-di NS 2 d i d j 

acts to the right. The current term is just 



(150) 



(151) 



(152) 



(153) 



S; 



(154) 



25 



where j M = (p, — j) . Hence 
where 



J [dA»] exp||z^M M ^-z^j^| 



-V z + a^ND 2 



(a^NS - I) D+Oi 



(a^NS - 1) D~di -nSij + (a^NS 2 - S) did. 



Assuming M„v has an inverse, we may integrate to find 



where 



The calculation of G goes as follows. We note ( |158|) is equivalent to 



[lvxVx— y^J y— z <J x—zi 



l.C 



J x 



Cp X M^ ux G x 5 . 



Assuming G v x has a DT fourier transform G v g , we may write 



e px M^ x <b e qx Gf = -<fj. 



Now e g:r is an eigenfunction of the operator M„ vx , so we may write 



where 



M, 



jivq 



q 2 + a~ 1 NuD 2 u i (1 - a^NuSu) D'uq 1 

% (1 - at^NuSu) D+uq* q 2 5 ij - (a- l NuS 2 u - Su) q l q j 

in block form, where q = (it, q). Hence we find 



M^ vp G v p x — —5*. 



Now if we write 



G 



a 


bp T ' 


b*p 


c + dpp T 


A 


B*p T 


Bp 


C + Dpp T 
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then a solution is 



.4 

C 



c + dp 2 
ac + p 2 A ' 



B 



D 



ac + p 2 A 
A 

c [ac + p 2 A] 



(166) 



where A = ad — b*b. Now from (|163|) we read off 



p 2 + a~ 1 NzD 2 z, b = i (1 - aT x NzSz) D' 1 ; 



P 



d = Sz — a 1 NzS 2 z 



(167) 



and now we may readily compute the components of G vX . 

The solution is complicated in the case of arbitrary N. We may now make a careful 
choice which reduces the complexity and brings us closer to the CT results. The particular 
choice which works here is 



which means that N (U n ) is the operator inverse of S (U n ) . 

Some concern may be expressed at this choice, as it involves step operators of all orders 
and may be undefined for some sequences. However, we can allay fears by noting that in the 
DT fourier transform space, the reciprocal of Sz = (z + 4 + 1/z) /6 exists, provided we are 
on the unit circle in the complex z plane. We may use this to define the inverse operator 
N (U n ) = S^ 1 (U n ) by its action in DT fourier transform space, viz 



s- 1 (U n ) f n 



z n Sz 



(169) 



We note that Sz has two zeros, located at z = —2 + v^3 and z = —2 — \^3. The first is inside 
the unit circle and contributes an additional simple pole away from the origin. Hence the 
above integral should be well defined and readily evaluated using the calculus of residues. 

With the above choice for Nz, the DT Landau gauge propagator is obtained by setting 
a = and then we find 



GT (P) 



p2p2 



-S 2 zp 2 
iSzD + zp 



iSzD zp T 



-p Sij — Szp l p J 

in block form. This tends to the correct CT Landau gauge propagator 



(170) 



in the CT limit 



GT (P) 



1 / .... pfp" 



P 



p 2 p 2 



-P 2 
~P°P 



T 
-p"p 



-p 2 5ij 



pip] 



0, Sz -> 1, 



-iD + z 



P°, 



-iD~z 



P°. 



(171) 



(172) 
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Here we have used the reparametrisation 



z = e 



(173) 



where 



6 = Tp° + (T 3 ) 



(174) 



as discussed in ||. 

The DT Feynman gauge propagator is obtained by setting a 
(Sz) 1 gives 



1 and choosing Nz 



G7 (p) = 



Sz T 
-5 



i.i 



which tends to the CT Feynman gauge propagator 



G7 (p) 



if 

p2 p2 



1 T 

-5. 



ij 



(175) 



(176) 



in the CT limit, as expected. 

Two comments on this last result are in order. 



1. Our theory has not assumed any metric structure to space-time. We note the appear- 
ance of a factor Sz in the — component in ( |175| ) . The precise relationship of our 
DT space-time and any DT analogue of the CT metric tensor in special and general 
relativity awaits further investigation; 

2. In the DT Feynman gauge we find the first functional derivatives with respect to the 
external sources of the vacuum functional gives 

(OouMO^j = ^-Z [j] = S nPn Z \j] , (177) 
which is consistent with the classical equation (|8~8D . 



8 Concluding remarks 

The stage is now set for the next paper in this series, which will be to develop the DT 
Feynman rules for QED and investigate the renormalization programme. One of the principal 
motivations of this series was the possibility that the introduction of a new scale T would 
allow a more convincing approach to the elimination of divergences than other formulations. 
Whether this works remains to be seen. We note that one of the early and successful 
regularisation methods, the Pauli-Villiars technique, employs the introduction of a very 
large mass scale. In some sense, this is the other side of the coin to our approach, which is 
to consider very small time scales. Either way, there is a suggestion or hint that something 
intimately involved with cosmological scales, either extremely large or extremely small, is at 
the core of the problems with field theory. 
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A A first order formulation 

In the first paper of this series H we interpreted Cadzow's equations of motion fllBD as an 
expression of momentum conservation at the node sites. The momentum in that context 
is clearly a node quantity. We shall call such a momentum a node momentum. In this 
section we consider the possibility of introducing link variables related to node variables 
via Legendre transformations and these we shall call link momenta. The node and link 
momentum concepts are different in DT mechanics, but it turns out that in the CT limit 
T —>■ they coincide. It is part of the problem of quantisation in DT mechanics that it 
is not clear a priori whether node or link momenta should be used in setting up canonical 
quantisation algebras such as the Weyl-Heisenberg algebra. In the previous papers of this 
series we avoided this question by using the DT analogue of the Schwinger action principle 
to obtain quantum commutation relations consistent with our DT mechanics. 

The reason we should consider link momenta is because there is more than one way to 
arrive at the Euler-Lagrange equations for configuration space variables in CT mechanics. 
The usual way is called the second order formulation, wherein we start with the Lagrangian 
L = L (q, q,t) and use Hamilton's principle directly to get the Euler-Lagrange equations 

jt (m) t It (178) 

Another route is via the so-called first order formulation, which is related to the Hamil- 
tonian phase space approach, but pretends to live in an extended configuration space. We 
define the conjugate momenta 

P 4 (179) 

construct the Hamiltonian H (p, q,t) , and then define the first order Lagrangian 

L = p ■ q-H. (180) 

Now we apply Hamilton's action principle to L and get the extended configuration space 
equations 

dL , . 

8p t < 182 > 

These equations are disguised versions of Hamilton's equations 
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. dH . dH . . 

q 7aF' P = "W (183) 

in phase space and are equivalent to the CT Euler-Lagrange equations obtained from the 
second order formulation. This approach to configuration space mechanics is the basis of 
the ADM formulation of quantum gravity [16|. We note that equations (|181 ) and (|182 ) 



are remarkably similar to the two sorts of DT equations of motion fl46|) and (|47| ) for node 
and link variables respectively This leads us to expect the DT analogue of phase space 
momentum to be a link variable, rather than a node variable. 

When it comes to gauge theories, the construction of link momentum for matter fields is 
altered in a subtle way related to the virtual paths used to maintain gauge invariance, but 
none of this subtlety is seen or required in the CT limit. 

For completeness, we show that we may rewrite the free electromagnetic equations in a 
first order formulation. We introduce the gauge invariant link variable iz n and define the 
system function 

Ft EE 7T n ■ D+A n - \<K n ■ 7T n + 7T„ • V0„ - \ (B nA • B nA ) . (184) 

Then the equations of motion are 

D-7z n + S n V xB„ = 0, (185) 

c 

V • 7r n = 0, (186) 

c 

w n = D+A n + V0 n , (187) 



c 



which are equivalent to the second order equations derived above. By inspection, we see that 
7r n = — E n ,which was to be expected from an analysis of the CT first order formulation. 

B Invariants of the motion 

Consider a symmetry of the system function, i.e. a transformation of the node and link 
variables 

AS^AZ + SAZ, fi^fi + Sfi (188) 
such that the system function is left unchanged. Then 

If = iAS (x) jJ^F" + S A° +1 (x, j^F" 

JV-1 r 




n=M 



Using the DT equations of motion (f46"l, 47) we find the invariant of the motion 



f f)T n f) T n 

T I d 3 *5A a n (x) - V • j- = C n+ \ (190) 
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which we refer to as a Maeda-Noether invariant |T7j . Although variations of the link variables 
are involved in the overall variation (and indeed are essential), partial derivatives with respect 
to link variables do not occur in the explicit construction of the above invariant. This 
underlines the fact that link variables are not equivalent to node variables in a dynamical 
sense. 
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